Extreme Nonlinear Optics in a Femtosecond Enhancement Cavity 
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Intrinsic to the process of high-order harmonic generation is the creation of plasma and the 
resulting spatiotemporal distortions of the driving laser pulse. Inside a high-finesse cavity where 
the driver pulse and gas medium are reused, this can lead to optical bistability of the cavity-plasma 
system, accumulated self-phase modulation of the intracavity pulse, and coupling to higher order 
cavity modes. We present an experimental and theoretical study of these effects and discuss their 
implications for power scaling of intracavity high-order harmonic generation and extreme ultraviolet 
frequency combs. 

PACS numbers: 42.65.Re, 32.20.Xx, 42.65. Pc, 42.65.Ky 



High-order harmonic generation (HHG) is initiated 
through field ionization of a gas medium, and thus high 
peak laser intensities are needed. This is typically re- 
alized with low repetition rate (< 100 kHz) amplified 
femtosecond laser systems producing high energy pulses 
(> 100 /xJ), with average powers up to tens of Watts. In 
an alternative method, the modes of a frequency comb 
and a high-finesse external cavity can be locked, lead- 
ing to > 100 MHz repetition rates and multi-kW average 
powers, while still supporting peak intensities sufficient 
for HHG [1, 2]. Efficient conversion of such high power 
to the extreme ultraviolet (XUV) holds great promise as 
a route to table top, high average brightness sources of 
coherent XUV light. Such a source could have a large 
scientific impact, particularly in applications where high 
repetition rate and/or temporal coherence is demanded, 
such as XUV frequency metrology Q and time resolved 
photoemission Q or photoionization coincidence [f| spec- 
troscopies. 

Since the first demonstrations in 2005 [2, Q, much 
progress has been made in the construction of high power 
frequency combs @, 0] and special optics for the efficient 
output-coupling of the generated HHG light from the 
cavity [H Q. Recently, intracavity powers as high as 18 
kW have been reported while maintaining femtosecond 
pulses jioj ]. However, while plasma refraction has long 
been known to be a critical issue in the phase matching 
of HHG, a cavity enhanced system operates with the ad- 
ditional constraints that the plasma must not ruin the 
cavity finesse or destroy the delicate resonance between 
comb and cavity. The effects of the plasma on the cavity 
performance have remained open questions and are the 
subject of this Letter. Only with these understandings 
have we been able to improve the cavity HHG yield by 
more than an order of magnitude and report here record 
levels of HHG performance in the XUV. Furthermore, 
understanding and controlling extreme-nonlinear-optical 
modulation of the intracavity light is critical for the phase 
coherence of the generated XUV frequency comb, be- 
cause phase modulation of the fundamental appears mul- 
tiplied by the harmonic order in the XUV light. Indeed, 
limited knowledge of the plasma phase shifts inherent to 



HHG was the dominant source of systematic error in a 
recent demonstration of XUV spectroscopy [U . 

The propagation of femtosecond pulses in a field- 
ionizing gas medium has been previously studied theoret- 
ically and experimentally lllll2j . The laser field loses en- 
ergy to the production of photoelectrons and experiences 
a rapid ramp in the index of refraction on the femtosec- 
ond time scale of the pulse, resulting in self-phase modu- 
lation and a frequency shift to the blue. Also, due to the 
high nonlinearity of field ionization, the plasma generated 
by a focused laser beam is far from spatially uniform, and 
acts as a negative lens, defocusing the beam. In a high- 
finesse cavity, the strength of these interactions between 
pulse and medium is effectively increased for two reasons. 
First, because the intracavity pulse is reused, pulse and 
wavefront distortion accumulates over many round trips. 
Second, if the repetition rate is sufficiently high, neither 
plasma expansion [13| nor gas jet flow is fast enough to 
clear the plasma from the focal volume and replenish it 
with neutral atoms within one cavity round trip. The in- 
tracavity pulse then gets an additional phase shift due to 
the residual plasma from previous round trips. In these 
ways, the plasma acts as a highly nonlinear optical ele- 
ment in the cavity, with responses on femtosecond and 
nanosecond time scales, presenting a novel and nontrivial 
problem in nonlinear optics. 

Our general cavity layout has been described previ- 
ously in Ref. @. The system is pumped with a Yb + 
fiber comb system capable of delivering up to 80 W of 
average power in 120 fs pulses with a repetition rate 
= 154 MHz and a center wavelength A c = 1070 nm 
Xe is introduced at the intracavity focus with a 200 
jxm diameter glass nozzle. The focus size (FWHM) was 
estimated to be 29 /xm (horizontal) x 17 /jm (vertical) 
from an ABCD matrix analysis of the cavity and inspec- 
tion of the frequency spacings of the cavity's higher order 
modes. The cavity is strongly overcoupled [T3|, with the 
round trip loss dominated by the 1.5% transmission of 
the input coupler, for a finesse of T ~ 400 and an op- 
timum power enhancement factor of ~ 250. In practice, 
a lower power enhancement of 200-220 is observed due 
to imperfect mode matching. The intracavity power is 
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FIG. 1: a) Intracavity power, b) PDH error signal, and c) 
13th harmonic power for a comb-cavity detuning scan with 
40-W pump power. The solid colored lines are experimental 
data and the black dashed lines are model results. The solid 
gray lines illustrate the intracavity power and error signal in 
the absence of gas. 



measured with a photodiode monitoring the transmitted 
light through a cavity mirror of calibrated transmission. 
Frequency-lock between the comb and cavity is accom- 
plished using the Pound-Drever-Hall (PDH) technique 
[15j | . The generated high order harmonics are coupled 
out of the cavity with a small period diffraction grating 
etched into the surface of a high reflector 0. The 13 th 
harmonic power is measured with a Si XUV photodiode 
coated with indium (IRD Inc. AXUV100In/MgF 2 ), and 
was confirmed to be insensitive to scattered IR and 3 rd 
harmonic light in the present geometry. 

Resonance between the pump comb light and the op- 
tical cavity is achieved when the optical frequencies of 
the comb lines are matched to the respective optical res- 
onances of the cavity. The nonlinear response of the sys- 
tem is clearly observed by sweeping the optical detun- 
ing Av = f com b — ^cavity across this resonance condition, 
as shown in Fig. [T] To acquire the data, the carrier- 
envelope offset frequency of the comb is first matched 
to that of cavity, and then the cavity length is scanned 
across the resonance while the intracavity power, PDH 
error signal, and 13th harmonic power are recorded with 
a digital oscilloscope. The detuning is scanned from neg- 
ative to positive Av in Fig. Q] and the frequency axis 
is calibrated using the resonances of the PDH sidebands 
observed in the error signal at ± 1.3 MHz. In the ab- 
sence of gas (gray curves), the intracavity power is an 
approximately Lorentzian function of detuning, and the 
PDH error signal, proportional to the first derivative of 
the line shape for small detuning, crosses zero at reso- 
nance. With gas flow and the onset of plasma creation, 



the intracavity power is reduced and the resonance con- 
dition is shifted to a blue detuning of more than a full 
cavity linewidth. 

All of the aforementioned consequences of light prop- 
agation in a field-ionizing medium can contribute to the 
nonlinear response of the cavity. To unravel their rel- 
ative contributions, we have developed a simple model 
appropriate for multicycle pulses under the present cir- 
cumstances, in which the cavity finesse is high and the 
gas medium density-length product sufficiently low that 
the change in the intracavity pulse induced in one round 
trip is small. We use a frequency domain approach to de- 
scribe the evolution of the intracavity comb on the time 
scale of tens of cavity round trips, denoted with variable 
T, and a time domain description to calculate the plasma 
ionization, loss, and high order harmonic generation on 
the femtosecond time scale, denoted with variable t. The 
expansion dynamics of laser produced Kr plasmas have 
been studied by Kanter et. 
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13j, who found the cen- 
terline density to decay with a time constant of ~ 20 ns 
governed by the ion speed of sound. The time scale for 
jet flow to replenish the focal volume is similar, on the or- 
der of 70 ns. The plasma decay is then fast compared to 
dynamics of the intracavity comb, with a characteristic 
time scale of F/(irf Tep ) ~ 800 ns, but slow on the fem- 
tosecond time scale of the pulse. Expressing the plasma 
density as por], where po is the atom number density and 
r] is the fraction of atoms that are ionized, the ionization 
fraction is then logically divided into two components: a 
steady-state component fj(T) that effectively tracks the 
intracavity pulse intensity, and a dynamic component 
r](t) responsible for self-phase modulation. 

The intracavity electric field is described by a discrete 
set of T-dependent spectral elements, Aj(T) = A(T,uij), 
each representing a portion of the intracavity spectrum 
with optical frequency Wj. The time evolution of the 
spectral elements is governed by a set of coupled differ- 
ential equations: 



dT 



+i(6 j +4> J )A d +f j ({A k }), (1) 



where 5 is the input coupler transmission, 7 is the round 
trip loss of the rest of the cavity, and Ej = E(ujj) is 
a set of spectral elements describing the electric field of 
the pump light. T is expressed in units of the cavity 
round trip time, l// rcp = 6.5 ns. The 0j describe the 
round trip phase shifts of the intracavity Aj with re- 
spect to the pump Ej due to the cavity detuning, anal- 
ogous to Av of Fig [1] The differential equations (1) 
are coupled through the steady-state plasma phase shifts 
(j)j(Ai, An) and the function fj(A\, ...,An), which ac- 
counts for the round-trip self-phase modulation and loss. 
In this way, different portions of the intracavity comb in- 
teract nonlinearly via the plasma. Both <f>j and fj are 
calculated from the t-domain description of the intracav- 
ity pulse as described below. 
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While the sub-cycle dynamics of field ionization, with 
very high frequency content, are essential for understand- 
ing HHG, they are not critical for explaining the distor- 
tion of the fundamental pulse if the pulse is multicycle 
and the ionization per cycle is small. Accordingly, at each 
time step, ionization, self phase modulation, and loss are 
calculated using an electric field envelope £(t), calculated 
from the Aj via discrete time Fourier transform (DFT). 
The empirical formula of Tong and Lin jl6| is used to 
construct cycle-averaged and peak ionization rates, w{£ ) 
and w{£). The ionization fraction during the pulse, r](t) 
is then calculated via 



?7(t) = l-[l-f?(T)]exp 



dt w(£) 



(2) 



The steady-state ionization fj(T) is calculated from a 
balance between the plasma creation per round trip, 
fj = rj(t = oo)—fj(T), and the simple decay law fj = —k p fj, 
where k p is a constant. 

To calculate the changes in the field envelope upon one 
pass through the plasma, we introduce the following ap- 
proximate expression derived from the propagation equa- 
tion of Geissler 111: 



A£(t) = ai£(t)\ -ir \ c p L[r,(t)-fj(T)} 



-I, 



2tt p Q L [1 - v (t)] w{£) 
\£\* 



(3) 



where L is the medium length, I p is the ionization po- 
tential of the target gas, ro is the classical radius of the 
electron, and c is the speed of light. The first term on 
the right-hand side accounts for the dynamic change in 
the plasma refraction during the pulse. Note that the 
steady state component, already included in <fij in equa- 
tion Eq. ([T]), has been subtracted off. The second term, 
essentially a low pass filtered version of the ionization 
loss term of Geissler 11| , accounts for the energy lost by 
the laser in ionizing the medium. Equation ([3]) neglects 
energy transferred to the electrons via their nonzero drift 
velocity, which is small compared to the ionization loss 
for Keldysh parameters of order unity or larger ll| . The 
coupling terms fj on the RHS of ([1]) are calculated from 
the inverse DFT of A£ . The reduction factor a.\ < 1 
crudely accounts for the 3D effect that only the center 
of the Gaussian focus experiences self-phase modulation, 
so that the effective modulation averaged over the entire 
wavefront is less. We include a similar factor oii < 1 in 
the calculation of 4>j from fj(T) via 

<t>j = -a 2 f]r Q X c Lp [l - Acjj/wc + (Acjj/w c ) 2 - ...] 

(4) 

where Aujj is the optical angular frequency difference 
from the carrier frequency u> c . For comparison with ex- 
periment, an HHG signal is estimated in the model by 
approximating the HHG dipole as d oc po(l — rj(t))y/w(£) 



and accounting for on-axis phase matching and absorp- 
tion as described in [17[. The model can also be easily 
extended to include the dynamic response of the system 
to perturbations with the frequency- lock servo loop en- 
gaged by adding additional differential equations to the 
coupled set JT]) [181 ] . 

With model parameters in the experimentally reason- 
able ranges a x , a 2 ~ 0.2 - 0.5, k p ~ 0.1 - 0.2, L = 200 
jttm, and po ~ 1 — 3 x 10 17 cm~ 3 , the model reproduces 
the measured data. The black dashed curves in Fig. Q] 
show model results from scanning the 6j from negative 
to positive, which simulates the cavity length scan. The 
model's error signal and HHG signal are scaled for com- 
parison with the experimental data, but the intracavity 
power is not. The observed distortion of the cavity reso- 
nance line shape can be understood as a new manifesta- 
tion of self- locking phenomena [2oT ] , where the non- linear 
phase shifts here are provided by the intracavity plasma. 
For large positive 9j, the systems exhibits optical bista- 
bility, with the coupled equations ([IJ having two pertur- 
bation stable steady-state solutions: a high intracavity 
power solution in which the plasma phase shifts, with 
the dominant contribution from <fij , compensate the cav- 
ity length detuning, and a low intracavity power solution 
with no gas ionization. In the frequency scan, self-locking 
is lost when the ionization of the medium saturates and 
the plasma phase shift cannot be increased further, caus- 
ing the system to jump to the low power solution. When 
actively locking the optical frequencies of the comb and 
cavity, high power operation is found to be perturbation 
unstable at the resonance condition (error signal = 0), 
and laser frequency or intensity noise can make the sys- 
tem jump to the low power solution. The frequency-lock 
servo then returns the system to resonance, producing os- 
cillations as shown in the inset of Fig. [21 with the servo 
loop effectively scanning through the self-locking curve 
of Fig. Q] repeatedly. These oscillations reduce the duty 
cycle and average power of HHG by a factor of ~ 2 and 
are deleterious to the coherence of the XUV comb. 

Stable locked operation can be maintained by intro- 
ducing an offset to the error signal such that the servo 
aims to maintain a frequency offset A between the shifted 
resonance condition of the cavity and the pumping comb. 
However, the intracavity pulse and spectrum are dis- 
torted from those of ideal linear operation. Fig. 2a) 
shows the measured intracavity spectra with and without 
target gas for A sa 40 kHz and 38 W pump power. The 
intracavity power is reduced by 11% and the spectrum 
is shifted to the blue. More severe problems with servo 
stability and reductions in the intracavity power occur at 
higher cavity finesse, reducing the attainable HHG power 
19]. Model results under similar conditions are shown in 
Figs. [2}d) and EJ:). In the time domain, the beginning of 
the intracavity pulse maintains better constructive inter- 
ference with the pump comb than the end of the pulse, 
which experiences a larger plasma phase shift. In the 
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FIG. 2: a) Measured intracavity spectra with (red) and with- 
out (gray) gas for offset A m 40 kHz and 38 W pump power. 
The intracavity power is reduced by 11%. b) The high power 
solution under similar conditions, c) Corresponding time do- 
main \£(t)\ and r)(t) (black dashed line, right scale). Inset: 
Intracavity power vs. time with (red curve) and without (blue 
curve) the lock offset A. With A = 0, the system oscillates 
on the ~ 20 (is time scales characteristic of the servo loop. 



frequency domain, the self-phase modulation appears as 
loss on the red side of the spectrum and gain for the blue, 
and the intracavity spectrum shows a corresponding shift 
to the blue. Repeating the calculation without the loss 
term in equation ([3]) yields nearly identical results, indi- 
cating that self phase modulation is the dominant effect 
limiting the intracavity power. 

While the ID model captures the essential physics, 
plasma lensing effects are present and manifest them- 
selves in the cavity-plasma system in a novel way. The 
spatial profile of the plasma has higher spatial frequency 
content than the fundamental TEMoo mode of the cavity 
and can thus couple power into the cavity's higher-order 
modes. This is clearly observed as modulation of the 
intracavity power at 21 MHz and its harmonics shown 
in Fig. [3] Inspection of the cavity's higher-order mode 
structure indicates the modulation comes from beating 
between the fundamental TEMoo mode of the cavity with 
the progression of higher order modes of even symmetry 
in the vertical direction. We hypothesize that mode cou- 
pling in the vertical direction is preferentially excited due 
to the smaller spot size and stronger plasma gradient in 
this dimension. The transition to multimode operation is 
observed to occur abruptly as the intracavity power and 



FIG. 3: RF spectra of the intracavity power signal with 1 MHz 
resolution bandwidth. The beat frequencies are labeled with 
the higher-order modes that have the corresponding frequency 
offset from the fundamental TEMoo- 



gas flow are increased. While the threshold conditions 
are typically beyond those found optimal for HHG, it is 
important to understand and control this effect, as such 
high frequency modulation severely impacts the phase co- 
herence of the XUV comb. The power coupled into higher 
order modes and the multi-mode appearance threshold 
depend dramatically on the cavity finesse, as shown in 
Fig. [3J where increasing the finesse by a factor of ~ 2.5 
reduced the appearance threshold from 7.6 kW to 5.3 kW 
and increased the power in higher order modes by more 
than 2 orders of magnitude. 

In conclusion, plasma optical bistability, cavity en- 
hanced spectral blueshifting, and multimode operation 
induced by plasma lensing must be examined in the de- 
sign of intracavity HHG systems. Stable locked operation 
can be maintained in the midst of a large steady-state in- 
tracavity plasma by careful servo loop design. Self-phase 
modulation remains as the primary obstacle to increasing 
the intracavity power, but modest drops in the finesse can 
allow large gains in obtainable intracavity power and sup- 
pression of multimode operation. Multimode operation 
can also be suppressed with an intracavity aperture. The 
optimization of HHG must also be considered, as seen in 
Fig. [TJ where the HHG signal goes through a maximum 
while the intracavity power continues to increase. With 
the current over-coupled cavity design with a finesse of 
~ 400, we have observed more than 20 itW out-coupled 
power levels under stable locked operation in the 13th 
harmonic. Assuming an outcoupling efficiency of 9% 
this corresponds to a generated power of 220 /xW, or a 
spectral brightness of roughly 4 x 10 16 photons/ (s mm 2 
mrad 2 (0.1%bandwidth)), approximately one tenth that 
of the XUV undulators of the Advanced Light Source 
Synchrotron at this wavelength [21] . 
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Yb + fiber laser. 

Note: Since Submission, we have become aware of re- 
lated work 22] 
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